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^ ■ ABSTRACT 
O ■ 

' We analyze longslit spectral observations of fourteen Wolf-Rayet galaxies from 

the sample of Schaerer, Contini & Pindao (1999). All 14 galaxies show broad 



o 



X 



Wolf-Rayet emission in the blue region of the spectrum, consisting of a blend 



^ ■ of NIIIA4640, CIIIA4650, CIVA4658, and HeIIA4686 emission lines, which is a 

spectral characteristic of WN stars. Broad CIVA5808 emission, termed the red 
bump, is detected in 9 galaxies and CIIIA5996 is detected in 6 galaxies. These 



^ . emission features are due to WC stars. We derive the numbers of late WN and 

early WC stars from the luminosity of the blue and red bumps, respectively. 
The number of O stars is estimated from the luminosity of the H/3 emission 
line, after subtracting the contribution of WR stars. The Schaerer & Vacca 
1998 (hereafter SV98) models predict that the number of WR stars relative 
to O stars, Nwr/Nq, increases with metallicity. For low metallicity galaxies, 
the results agree with predictions of evolutionary synthesis models for galaxies 
with a burst of star formation, and indicates an IMF slope -2 < F < — 2.35 
in the low metallicity regime. For high metallicity galaxies our observations 
suggest a Salpeter IMF (F = —2.35) and an extended short burst. The main 
possible sources of error are the adopted luminosities for single WCE and 
WNL stars. We also report, for the first time, NGC 450 as a galaxy with 
WR characteristics. For NGC 450, we estimate the number of WN and WC 
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stars. The number ratio Nwr/Nq, and the equivalent widths of the blue 
bump, EW\4Qse, a-iid of the red bump, EWx^sos in NGC 450 are also in good 
agreement with the instantaneous burst model prediction for WR galaxies. 

Key words: galaxies: Starburst — galaxies: abundances — galaxies: evolu- 
tion. 

1 INTRODUCTION 

Wolf-Rayet galaxies are extragalactic objects whose spectra show direct signatures similar to 
those observed in Wolf-Rayet (WR) stars. The most common characteristic is the presence 
of a broad HeIIA4686 feature (the blue bump) originating in the stellar winds of WR stars 
(Schaerer Contini & Pindao 1999, hereafter SCP99). WR galaxies have long been known, 
with the first discovery of such spectral features in the blue compact galaxy He 2-10 (Allen, 
Wright & Goss 1976). The concept of WR galaxies was introduced by Osterbrock & Cohen 
(1982, hereafter OC82) and Conti (1991). 

The blue WR bump is often blended with nearby nebular emission lines of He, Fe, or Ar, 
and can show several broad stellar emission components (NHIA4640, CHIA4650, HeHA4686) 
which are difficult to deblend in most low- or medium-resolution spectra. These features 
originate in WR stars of WN and/or WC subtypes (OC82 and Conti 1991). The strongest 
emission line in WC stars is C1VA5808, which is very weak in WN stars. This "red W-R 
bump" has only rarely been observed. Where the data is available, CIVA5808 is generally 
weaker than HeIIA4686. 

In more distant galaxies, WR stars can only be indirectly detected, by observing the 
integrated spectra of the galaxies. Strong star formation activity indicates the presence of 
a large number of massive stars, most of which evolve through the WR phase. At a given 
stage of the starburst, many WR stars appear, but only for a brief duration. Thus, the 
presence of WR features in these galaxies indicates recent star formation (<10 Myr) as well 
as the presence of massive stars {Minuiai > 25Mq) (Schaerer et al. 1999). This provides 
interesting constraints on recent star formation episodes in these objects (Maeder & Conti 
1994). Furthermore, metallicity plays an important role in regulating the lower mass limit 
above which a star passes through the WR phase. 

Despite their small number compared to other massive stars, especially in low-mctallicity 
galaxies, WR stars are numerous enough for their integrated emission to be detected. In 
this paper, we shall follow OC82 and Conti (1991): A WR galaxy is classified as such if 
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its integrated spectrum shows detectable WR broad features emitted by umesolved stellar 
clusters. 

The compilation of Conti (1991) included only 37 objects. Since then, the number of 
cataloged WR galaxies has increased rapidly, with more than 130 known today (Guseva, 
Izotov & Thuan 2000, hereafter GITOO; SCP99). WR galaxies do not form a homogeneous 
class, exhibiting a variety of morphologies. Among the WR galaxies we find low-mass blue 
compact dwarf (BCDs), irregular galaxies, massive spirals and ultra-luminous merging IRAS 
galaxies. Recent studies show that WR features are also seen in LINERs and Seyfert 2 
galaxies (OC82; Ho et al. 1995; Heckman et al. 1997; Schmitt et al. 1998 and Contini et al. 
2001). The possibility of detecting WR stars in central cluster galaxies out to a redshift of 
z > 0.25 is discussed in Allen (1995). 

The number of WR stars relative to massive stars is highly dependent on metallicity. 
Theoretical evolutionary models predict that at fixed metallicity, the ratio between WR and 
other massive stars varies strongly with the age of the starburst (Mass-Hesse & Kunth 1991; 
Maeder 1991; Maeder & Meynet 1994; Meynet 1995; SV98). The maximum value of this 
ratio decreases from 1 to 0.02 when the metallicity decreases from Zq to Zq/50 (GITOO). 
Similarly, the duration of the WR stage in the starburst also decreases with decreasing 
metallicity. Hence, the number of galaxies with extremely low metallicity containing WR 
stellar populations is expected to be small. 

GITOO derived the number of WCE and WNL stars from the luminosity of the red and 
blue bumps, respectively, and the number of O stars from the H/? luminosity, for 39 WR 
galaxies with heavy element mass fractions between 1/50 and twice solar. In their sample, 
the blue bump consists of an unresolved blend of WR and nebular lines. They proposed a new 
technique to derive the number of WNL stars using NIIIA4512 and SiIIIA4565 emission hues. 
They found that the relative number of Wolf-Rayet stars ^wRf^io+WR) and Nwc/^wn 
derived from observations arc in satisfactory agreement with theoretical predictions (SV98). 
The results obtained for the extremely metal-poor galaxies disagree with model predictions 
likely due to the low emission line luminosity for WCE stars in metal-poor models. 

Schaerer et al. (2000, hereafter SGITOO), using five metal rich objects from GITOO and 
new results on Mrk 309, attempted to constrain the properties of massive star populations 
and star formation histories by comparing their observations with evolutionary synthesis 
models. They found that extended burst durations of ~ 4 — 10 Myr or a superposition of 
several bursts were required to produce the observed WR population and red supergiant 
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features. The burst durations are longer than those obtained for other objects in Schaerer 
et aL (1999) using the same models. 

Pindao et al. (2002) analyzed the spectra of 85 high metallicity disk HII regions of nearby 
spiral galaxies. In contrast with previous studies of low metallicity galaxies, they found 
smaller values of I(WR)/I(H/3) than predicted from evolutionary models at corresponding 
metallicities. They suggested the use of two WR luminosity regimes to correct the model 
predictions. 

The goal of our study of these 14 galaxies is to search for and confirm the presence of WN 
and WC stars in galaxies with different metallicities, and to compare the results obtained 
for this sample with predictions from evolutionary synthesis models (SV98) and starbust99 
(Leitherer et al. 1999). We also report NGC 450 as a newly classified WR galaxy. 

The paper is structured as follows: the observations and procedures used to reduce the 
data are described in §2. In this section we also discuss corrections for reddening and un- 
derlying absorption affecting the emission lines. Section 3 describes how contamination of 
a starburst spectrum by the presence of SN He and AGN is accounted for, while Section 4 
explains how the physical parameters of the gas are calculated. The massive star population 
and constraints on the evolutionary tracks of the starburst regions (age, burst duration, and 
IMF) are derived in §5 from a comparison with evolutionary synthesis models. Finally, our 
main results are summarized and discussed in §6. 

2 SPECTROSCOPIC OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

We observed 14 galaxies from the sample of Wolf-Rayet galaxies and extragalactic HII regions 
presented in SCP99. NGC 6764 was observed at two position angles, 67° and 90°, with the 
former along the major axis of the galaxy and the latter spanning a secondary emission 
region near the central region. We suspected that WR stars may be present in this region, 
but the spectra do not confirm this hypothesis. 

Data were gathered at two different sites. We used the Palomar 200-inch telescope with 
the longslit Double Spectrograph (Oke & Gunn 1982) on UT 1999 October 10-11, with the 
1200 1/mm grating, blazed at 5000A, yielding a pixel size of 0.62" x I.SA in the blue and 
0.47" X I.7A in the red, and a total wavelength coverage of 3600 - 6700A. The slit width 
was set to 1", resulting in a spectral resolution of ~ 5.6A in the blue and ~ 5.7A in the red. 
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The slit length was 180". Additional data were taken at the 3.6m ESO/NTT telescope using 
the ESO Multi-Mode Instrument (EMMI) with CCD36 on UT 1999 April 17. The total 
wavelength coverage was 4000-6600A. The slit width was set to 1", resulting in a spectral 
resolution of ~ 5.9A. The length of the slit was 120". 

The sht was ahgned along the major axis of the galaxy when possible, and centered on 
the brightest region of the target. The journal of observations of all objects in the sample 
is provided in Table 1. Figures l(a)-(o) show the images of each galaxy taken from the DSS 
with the slit position overlaid. 

2.2 Data Reduction 

Initial data reduction was carried out following standard procedures with the IRAF task 
ccdproc. The package includes bias subtraction, fiat-field correction, subtraction of the night 
sky background and bad column removal. The bias level was subtracted from each frame 
using the overscan region of the CCD chip. 

Spectrum extraction was performed using the IRAF task longslit. The task includes 
correction for atmospheric extinction, wavelength calibration and flux calibration. Cosmic 
rays were removed using the cosmicrays task with a threshold of 5%. 

Spectrophotometric standard stars were observed each night to perform flux calibration. 
Arc lamps were taken before and after each exposure in order to provide accurate wavelength 
calibration. An average rms of 0.1 A was obtained for the pixel to wavelength fit using a 
third order spline. 

We use a standard extraction aperture for each object whose width is such that the peak 
intensity of the H/3 line decreases by 80% from the center, along the spatial direction. This 

width is set independently for each side of the peak. For objects with more than one emission 
region, different apertures were set for each knot. With this size we obtain the necessary 
spectral properties for this study while minimizing possible contamination from the adjacent 
stellar population. 

The signal to noise ratio, S/N, was determined assuming Poisson statistics and using the 
readout noise and gain of the CCD, the number of combined spectra, and the sky value. For 
NTT data the readout noise is 5.43e~^ with a gain of 2.18 e~^ per ADU. For the Palomar 
observations, the readout noise is 8.6 e^^ with a gain of 2.13 e~^ per ADU in the blue channel 
and 7.5 e"^ and 2.00 e"^ per ADU in the red channel. The values of S/N are given in Table 
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Table 1. Journal of Observations 



Galaxy 


a 


5 


Exp. Time 


S/N 


P.A. Observat. Type 


D 


Scale mag 




^hh .mm .ss ^ 


(°:':") 


(s) 




n 


(Mpc) 


(pc/") (Band) 



UM 48 


0:36:10 


4:38:06 


3x 


1200 


96/73 


43 


Palomar 


S 


65.22 


316 


13.35(J) 


NGC 450 


1:15:31 


-0:51:38 


1500 


+ 1200 


90/77 


80 


Palomar 


SAB(s)cd 


22.46 


109 


11.40(J) 


MRK 712 


9:56:42 


15:38:04 


4x 


900 


84 


108 


NTT 


SBbc 


55.87 


271 


14.47(B) 


MRK 1271 


10:56:08 


6:10:15 


4x 


900 


80 


98 


NTT 


Compact 


8.79 


43 


14.80(mP) 


NGC 4385 


12:25:42 


0:34:23 


4x 


900 


96 


192 


NTT 


SB{rs)0+ 


23.35 


113 


10.64(J) 


NGC 4861 


12:59:00 


34:50:43 


1800 


+ 600 


160/97 


19 


Palomar 


SB(s)m 


10.93 


53 


12.44(J) 


NGC 5430 


14:00:46 


59:19:35 


2x 


1200 


140/104 


146 


Palomar 


SB(s)b 


40.55 


197 


lO.Ol(J) 


NGC 5471 


14:04:29 


54:23:49 


2x 


1200 


91/75 


127 


Palomar 


HII 


3.86 


19 


14.76(J) 


MRK 475 


14:39:05 


23:37:28 


2x 


1800 


84/70 


80 


Palomar 


BCD 


7.96 


39 


15.46(B) 


Fairall 44 


18:13:38 


-57:43:59 


2x 


900 


90 


105 


NTT 


S? pec 


63.62 


308 


12.34(J) 


NGC 6764 


19:08:15 


50:55:57 


2x 


1200 


125/95 


67/90 


Palomar 


SB(s)bc 


31.26 


152 


10.56(J) 


MRK 309 


22:52:34 


24:43:50 


2x 1500 + 1800 


123/90 


140 


Palomar 


Sa 


167.09 


810 


13.05(J) 


III ZW 107 


23:30:09 


25:32:02 


1800 


+ 1200 


125/95 


175 


Palomar 


Sb 


76.50 


371 


14.06(J) 


NGC 7714 


23:36:13 


2:09:21 


2x 


1200 


99/71 


122 


Palomar 


SB(s)b:pec 


37.05 


180 


10.77(J) 



(a) and (5) coordinates (,12000) arc for the slit center. The distance (D) is obtained from the redshift measured directly from 
the object spectrum using Ho =75 km Mpc^^. NTT and Palomar indicate the observatory where the object was observed. 
The object type is obtained from NED. The linear scale is obtained using the computed distance, D. *For more details on the 
quoted magnitudes we refer the reader to the NED-IPAC database. 

1 for the NTT, as well as the blue and red parts of the Palomar spectra, where they are 
averages over a range of 100 A around 4686A and 5878A, respectively. 

2.2.1 Reddening Correction and Underlying Balmer Absorption 

The reddening correction is obtained from the Balmer line ratios, using the extinction law 
of Cardelli, Clayton & Mathis (1989) (hereafter CCM89), and the theoretical Balmer emis- 
sion line ratios for case B recombination: I{Ha)/I{Hj3) = 2.88 and I{H'~f)/I{Hf3) = 0.47 
(Brocklehurst 1971). We obtain the color excess, E{B — V), using the theoretical ratio of 
Balmer to F[Ha)/F[Hj3) measurements, according to 

(F{Ha)\ 



E{B - y) = 2.33 X 



(1) 



In three cases (Mrk 309, Mrk 1271, and Fairall 44) for which F{Ha) is not available, we 
use solely the measured F{H^) / F{Hl3) ratio to obtain the reddening correction using 



E{B -V)= 4.71 X 



(2) 



The resulting color excess is used to correct each spectrum with the IRAF task deredden. 
E(B-V) is easily converted to cHP using the relation E{B -V) 0.677 x cHP (Vogel et al. 
1993). 

In Tabic 2 wc compare the reddening E(B-V) measured in this work with that obtained 
by other authors. Furthermore, wc have compared our derived color excesses with other 
measurements taken from the literature in Figure 2. The top panel shows this comparison, 
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Table 2. E(B-V): Comparison with Literature. 





()l>scrv(-(l 


Lilci al lire 


l^(!r('l(!liCC~ 


Fairall44 


0.47 


0.60 


5 


III Zw 107 


0.53 


0.47 


9 


Mrk309 


0.73 


0.71, 0.63, 0.70 


9, 2, 12 


Mrk475 


0.21 


0.10 


7 


Mrk712 


0.03 


0.18 


3 


A/Tt-VI 971 


u. ±o 


n in 


Q 

O 


NGC450 


0.55 






NGC4385 


0.64 


0.60, 0.52 


2, 14 


NGC4861 


0.43 


0.20, (0.49, 0.62, 0.77) 


2,(1) 


NGC5430 


0.67 


1.14 


4 


NGC5471 


0.21 


0.15 


13 


NGC6764 


0.82,0.80 


0.61,(0.87) 


10, (6) 


NGC7714 


0.40 


(0.39,041), 0.43 


(2), 11 


UM48 


0.57 


0.49 


15 



(1) Barth et al. 1994, (2) Calzetti 1997, (3) Contini ct al. 1995, (4) Contini ct al. 1997, (5) Durrct 1990, (6) Eckart ct al. 1996, 
(7) Izotov et al. 1994, (8) Izotov & Thuan 1998, (9) Gil de Paz et al. 2000, (10) Keel 1982, (11) Mattila & Meikle 2001, (12) 
Osterbrock & Cohen 1982, (13) Rosa & Bevenuti 1994, (14) Salzer et al. 1989 ,(15) Vogel et al. 1993 

while the bottom panel shows the deviations as a function of E(B-V). The median difference 
is 0.06 with a standard deviation of 0.18. 

We note some differences between the values obtained for some of the objects. The 
spectroscopic results taken from the literature are inhomogeneous in their instrumental 
accuracy, observational techniques, and the S/N ratios of individual observations. However, it 
is unlikely that instrumental differences dominate over variations in measured lines intensities 
due to observing emission from different locations within each galaxy or star formation 
region, different sized apertures, or slightly different position angles used by the authors. 

After correcting the spectra for the total interstellar extinction using the Balmer line 
ratios, some objects show F{H^)/F{HP) ratios lower than the theoretically predicted val- 
ues. This effect is more significant for objects where the H/3 equivalent width is less than 
70A (McCall, Rybski & Shields 1985, hereafter MRS85). The consequence of such under- 
lying absorption is an overestimate of the reddening correction, which is then reflected in 
all the observed emission lines. To correct for this effect, we use the method proposed by 
MRS85 to compute the equivalent width in absorption, EWabs, given by the expression, 

( F{Hj) \ ^ I{Hj){l + EWaks/EW{HP)Y-' 
I{HP) il + EWa,jEW{Hj)) " 

The parameter e is defined by 

(AiHjX _ AiH0) \ 
^ ( A(Ha) _ A(Hff) \ 

\ Ay Ay ) 

and represents the extinction reddening law (CCM89 and MRS85). 
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In the limit EWabs = 0, iF{H'j) / F{H I3))dered reaches the theoretical ratio (/(i?7) /I (Hp) = 
0.47). In the opposite limit, as EWabs approaches unity, the ratio {F{H'y)/F{HP))dered is 
significantly altered from the theoretical expectation. The corrected equivalent widths of the 
Balmer fines are obtained by adding EWabs to the observed equivalent widths. 

The correction in the extinction due to underlying absorption lines, E{B-Y) abs, is given 
by the expression 

1 + {EWabs/ EW{Ha)) \ . . 

l + {EWabs/EW{Hp))) ^ ^ 

The value of EWabs for each object is listed in Table 3, along with the equivalent widths of 

the emission lines EW{Ha), EW{HP) and EW{Hj). Reddening and underlying absorption 

corrected intensities are also given in the same Table. 



E{B -V)^ -2.33 X Log 



3 HARD IONIZATION SOURCES 

The contamination of a starburst by nearby objects such as SN He remnants and AGN 
could occur, producing Wolf-Rayet-like features in the spectrum (Masegosa et al. 1991). 
Embedded or superimposed supernova remnants (SNRs) can strongly affect the emission- 
line properties in a predominantly photoionized HII region {e.g., Peimbert, Sarmiento, & 
Fierro 1991). It is crucial to distinguish between hot stars and other ionizing photon sources. 
Prom an observational viewpoint, the presence of shocks in an ionized nebula can be detected 
by the enhancement of low excitation lines, in particular [01], [SII] and [Nil] (Masegosa et 
ah 1991). 

We used Veilleux & Osterbrock (1987) diagnostic diagrams to verify the ionization 

sources in our galaxy sample. Since [SII]AA6717, 6731 are not available in our spectra, we 

use diagnostic diagrams plotting [OI]A6300/Hq; and [NII]A6584/Hq; fine ratios against the 

[OIII]A5007/H/3 ratios (Pigure 3 and 4). 

The theoretical starburst limits shown in Pigure 3 and 4 can be parameterized as: 

/og([Q///]A5007) ^ OM 

HP log{[NII]X658A/Ha) - 0.47 ^ ^ 

log{[OIII]Xm7 ^ 0J3 ^ 

HP log{[OI]\6300/Ha) + 0.59 ' ^ ' 

given by Kewley et al. (2001). The error range of their modeling in both planes is ± 0.1 dex. 

All objects in our sample arc located in the starburst and HII region loci in the plane 

of [OI]A6300/Ha and [NII]A6584/Ha versus [OIII]A5007/H/3 (Veilleux & Osterbrock 1987), 

suggesting that the main sources of ionizing photons are stellar, and not AGN. 

© 2004 RAS, MNRAS 000, 1-25 



Massive Star Populations in Wolf-Rayet Galaxies 9 



Table 3. Corrected Emission Line Intensities Relative to 113 = 1000 





Fairall44 


IIIZW107 


IVIRK309 


MRK475 


Mrk712 


Mrkl271 


[OII]A3727 


- 


2213±12 


703±98 


1093±4 


- 


- 


[NcIII]A3869 


- 


259±5 


- 


362.9±2.7 


- 


- 


H7 


475±2 


470±5 


466±46 


470±3 


470±3 


470±8 


[OIII]A4363 


26±9 


13±4 


- 


85±3 


14±2 


61±12 


HcIA4471 


53±9 


49±5 


59±38 


42±2 


67±2 


61±15 


[FeIII]A4656 


29±5 


13±7 


131±29 


5±3 


14±3 


26±7 


HeII4686 


19±5 


3±9 


- 


17±4 


- 


35±11 


[ArIV]A4711 


- 


14±10 


- 


8±3 


8±2 


- 


[ArlV]A4740 


- 


- 


- 


7±3 


- 


- 


HcIA4922 


- 


12it8 


- 


11±3 


19±2 


- 


[OIII]A4959 


710±5 


1215±7 


104±25 


1781±5 


1613±4 


1115±12 


[OIII]A5007 


2061±9 


3632±18 


381±33 


5326±15 


4881±11 


3096±29 


[NI]A5199 


- 


31±7 


- 


4±3 


12±3 


- 


[FeIII]A5271 


- 


- 


- 


7±3 


10±3 


- 


[C1III]A5518 


- 


10it2 


- 


8it2 


8±5 


- 






3±1 




4±2 


6it5 




[NII]A5754 




9±2 


38±25 


5±3 


7±3 




HeIA5876 


126±4 


129±2 


24±17 


99±1 


228±3 


119±9 


[OI]A6300 


48±6 


53±2 


26±22 


24±1 


36±3 


65±11 


[SIII]A6313 


8±4 


11±2 


_ 


24±1 


24±3 


15±10 


[OI]A6364 


18±4 


18±3 


_ 


9±1 


12±3 


72±11 


[NII]A6548 


_ 


98±2 


_ 


33±2 


268±2 


84±11 


Hq 


- 


2855±14 


- 


2870±8 


- 


2899±25 


[NII]A6584 




280±2 




66±1 




221±9 


EW(H7)(A) 


32.17±0.23 


26.65±0.35 


3.00±0.18 


73.08±3.35 


45.98±0.21 


15.42±0.20 


EW(H/3){A) 


67.64±0.26 


65.48±0.46 


10.18±0.79 


140.50±0.37 


139.10±0.29 


33.11±0.28 


EW(Ha)(l) 




329.60 ±3.70 




936.20±10.98 




109.20±1.16 


EW„fc,(A) 




-0.98 




-1.18 




-0.69 


F{H/3)(10-14 erg s"! cm-2) 


3.58±0.01 


12.77±0.06 


12.8±1.04 


7.52±0.03 


2.45±0.02 


0.72±0.01 


E(B-V) 


0.47 


0.53 


0.73 


0.21 


0.03 


0.15 


E(B-V)„6, 




-0.012 




-0.007 




-0.015 


Broad Emission Lines 


NIIIA4640 




7±5 


167±33 




82±10 




HeIIA4886 


59±9 


61±26 


115±20 


69±16 


98±7 


103±41 


CIIIA5696 






15±4 








CIVA5808 








36±7 






EWA4640(A) 




0.37±0.38 


1.36±0.21 




7.23±0.98 




EWA4686(A) 


3.47±0.56 


3.60±2.86 


0.97±0.20 


9.51±2.35 


8.67±0.74 


3.06±1.21 


EWA5696(A) 






0.26±0.05 








EWA5808(A) 








8.02±1.46 







4 GAS CHEMICAL ABUNDANCES 



The oxygen abundance 0/H is a crucial parameter for our study since one of the main 
purposes of this paper is to study the massive star populations in galaxies with different 
metallicities. In order to obtain the chemical abundances of the galaxies in our sample, 
three subclasses were defined depending on whether or not the [OIII]A4363, [OII]A3727 and 
[OIII]AA4959, 5007 features are present (see Kobulnicky et al. 1999): 

Class 1. [OIII]A4363, [OII]A3727 and [OIII]AA4959, 5007 are detected. When [OIII]A4363 is 
detected, the temperature is obtained from the [OIII]AA5007, 4959/A4363 line ratio, while 
the density is measured from the [SII]A6731/[SII]A6717 line ratio (McCall 1984). Since our 
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Table 3. Continued 





NGC450 


NGC4385 


NGC4861 


NGC5430 


NGC5471 


NGC6764 




1615±3 




626±2 




919±2 


1933it29 


[NcliI]A3869 


240±2 


- 


408±2 




399±2 


68±18 


H7 


481±2 


471±3 


470±2 


470ib23 


470±2 


470±12 


[OIII]A4363 


27±2 


13±3 


89±2 




66±2 


21±8 


HeIA4471 


48±2 


30±6 


37±2 


46±9 


39±2 


47±13 


[Fclll] A465d 


9±3 


35±6 


8±3 


49ib23 


26±4 


69±22 


HcIIA4686 


~ 


14±4 


8±2 


32±10 


25±1 


54±15 


[ArIVlA4711 


5±2 


- 


11±2 




14±1 


- 


[ArIVJA4740 






7±2 








HeIA4922 


13±2 


- 


5±2 


- 


9±1 


- 


[UiliJ A4959 


13d8±3 


231±3 


1817±4 


129±10 


1850±30 


168it6 


[Oin]A5007 


3864±7 


678±4 


5521±11 


375±12 


5583±90 


484±8 


[NIJA5199 


9±2 


19±3 


3±3 


12±12 


- 


88±14 


[FcIII]A5271 


6±3 


- 


- 




2±2 


- 


[C1III]A5518 


6±1 




4±2 


■ 


5±1 


7±5 


[C1III]A5538 


5±1 




3±2 


2±2 


6±2 




[NII1A5754 


2±1 


10±3 




12±4 


3±2 


29±7 


HeIA5876 


131±1 


127±3 


92±1 


119±5 


lllitl 


123±6 


[OI]A6300 


25±1 


25±4 


16±2 


- 


22±1 


159±7 


[SIII]A6313 


15±1 


10±4 


17±2 


6.7±3.3 


18.6±0.7 


- 


[OI]A6364 


9±1 


6±5 


5±1 


4±3 


8±1 


67±8 


[NII]A6548 


71±1 


503±4 


23±2 


441±11 


46±1 


726±12 


Hq 


2872±5 


2884±8 


28.5.5±22 


2837±49 


2866±5 


2806±35 


[iNilJAOOo4 


10 ( ±1 


14oy±0 


A V-LO 
4 I ±Z 




on 1 -1 




EW(H7)(A) 


93.61±0.31 


13.85±0.09 


48.25±0.32 


13.21±0.45 


51.68±0.18 


6.58±0.14 


EW(H/3){A) 


320.80±0.59 


42.78±0.11 


156.80±0.29 


40.85±0.41 


181.20it0.29 


22.46±0.27 


EW(Ha){A) 


2928.00±10.59 


199.10±1.66 


938.00±5.61 


264.60±1.25 


1028.00±6.35 


136.40±0.90 


EW„6,(A) 


-0.27 


-2.18 


0.00 


-0.05 


-0.40 


-0.72 


F(H/3)(10-" erg s"! cm-2) 


30.78±0.06 


48.63±0.15 


103.60±0.20 


37.17±0.63 


7.73±0.01 


59.81it0.83 


E(B-V) 


0.55 


0.64 


0.43 


0.67 


0.21 


0.82 


E(B-V)„i„ 


-0.001 


-0.039 


0.000 


-0.001 


-0.002 


-0.027 


Broad Emission Lines 


NIIIA4640 


21±16 


95±11 




127±26 




161±33 


HeIIA4886 


31±14 


103±8 


55±13 


105±18 


22±4 


63±21 


CIIIA5696 




31±31 


3±2 






10±5 


CIVA5808 


23±3 


62±34 


11±6 


41±10 




74±27 


EWA4640(A) 


5.14±3.80 


3.53±0.39 




5.50±0.94 




2.94±0.61 


EWA4686(A) 


7.61±3.80 


3.91±0.26 


7.08±1.75 


3.61±0.78 


3.14±0.69 


1.15±0.57 


EWA5696(A) 




2.16±2.12 


0.81±0.47 






0.30±0.18 


EWA5808(A) 


10.58±1.60 


4.55±2.41 


3.00±1.49 


3.16±1.77 




3.90±0.76 



spectra do not include the [SII] AA6717, 6731 lines, we are unable to measure the electron 
density, and we assume a fiducial value of 150 cm~^ for ng. 

The accuracy of this empirical method depends strongly on an accurate estimate of the 
electron temperature of the gas (Steigman, Viegas & Gruenwald 1997). The uncertainty in 
the electron temperature of the gas is directly related to the measurement of the [OIII]A4363 
line intensity. The oxygen abundance can be obtained by the empirical method first proposed 
by Peimbert & Costero (1969), which requires knowledge of the gas temperature and density. 
Class 2. [0I11]A4363 is not detected, while [011]A3727 and [0111]AA4959, 5007 arc detected. 
For those galaxies where the [OIII]A4363 auroral line is not detected, the oxygen abundance 
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Table 3. Continued 





NGC6764P.A.90 


NGC7714A 


NGC7714B 


UM48 


[Oil] \6lll 


17U3±31 


1 rtTr\ 1 

1970±o 


rii nT 1 1 n 

2107±10 


2152±29 


[Ne III]A3869 


130±23 


95±3 


280±4 


lllitlS 


H7 


471±15 


472.4±3.7 


469.9±4.4 


472±16 


[0 III]A4363 


19±8 


13.6±4.2 


39.3±4 


12±6 


HeIA4471 


59.7±16.5 


39.6±3.5 


46±4 


54±18 


[Fclll] A4656 


85±26 


30±5 


28±9 


35±23 


TT TT \ Anon 


58±23 


6±1 


4±4 


12±16 


r A 7"! \ ATI 1 

[AriVJA4711 






22±7 




[ArlVJ474U 










HeIA4922 




- 


12±4 


- 


r/~^TTTl \ Anr^f\ 

[OIII] A4959 


171±8 


496it40 


1180±6 


518±14 


Fi^TTTl \ 1- r\r\T 

[UillJA5U07 


484±11 


155U±140 


3529±14 


15o3±20 


[NI]A5199 


77±15 


18±4 


- 


- 


[FcIII]A5271 




27±4 


23±4 


- 


[C1III]A5518 


10.7±9 


7±4 






[C1III]A5538 


4.2±6 


7±4 






[NII1A5754 


20±7 


11±3 


llitlO 


- 


HeIA5876 


99±7 


125±3 


117±3 


125±5 


[OI]A6300 


156±8 


44±4 


46±4 


34±10 


[SIII]A6313 


- 


8±4 


15±4 


15±11 


[OI]A6364 


58±8 


14±6 


16±4 


13±7 


[NII]A6548 


708±13 


362±4 


131±4 


197±5 


Ho 


2811±40 


2823±123 


2872.2±12 


2869±31 


[iNilJAD0o4 


1 ni 1 OQ 
ly Izztzo 


OT"> 1 r 


QQ /I -1-/1 




EW(H7)(A) 


6.39±0.17 


12.13±0.08 


68.88±0.57 


9.95±0.31 


EW(H/3){A) 


22.04±0.30 


32.99±0.10 


191.30±0.72 


29.25±0.31 


EW(Ha){A) 


137.70±0.97 


98.48±0.63 


12.00±0.01 


156.30±1.09 


EW,b,(A) 


-0.94 


-1.44 


-0.18 


-1.29 


F(H/3)(10-" erg s"! cm-2) 


48.98±0.80 


86.15±0.32 


7.79±0.03 


9.32±0.12 


E(B-V) 


0.80 


0.52 


0.40 


0.57 


E(B-V)a(„ 


-0.036 


-0.029 


-0.014 


-0.035 


Broad Emission Lines 


NIIIA4640 


161±39 


41±9 






HeIIA4886 


64±23 


22±10 


54±14 


56±33 


CIIIA5696 


10.±7 


5.0±0.6 






CIVA5808 


81±29 


34±14 




22±12 


EWA4640(A) 


2.94±0.71 


1.41±0.25 






EWA4686(A) 


1.20±0.54 


0.60±0.27 


10.34±2.92 


1.50±0.84 


EWA5696(A) 


0.23±0.17 


0.22±0.03 






EWA5808(A) 


3.65±0.69 


1.77±0.08 




1.52±0.71 



is obtained from the empirical calibration suggested by McGaugh (1991) using the known 
R23 ratio, ([OII]A3727+[OIII]AA5007,4959)/iJ/3. 

The relation between oxygen abundance and R23 is double valued, but this degeneracy 
can be broken. We use the relation between [NII]A6584/[OII]A3727 and [NII]A6584/Hq; ratios 
to break the 0/H vs. R23 degeneracy (Contini et al. 2002). 

Kobulnicky et al. (1999) provide a polynomial fit to both metal-poor (lower branch), 

12 + log{0/H)i = 12 - 4.944 + 0.767a; + 0.0602a;2 - 7/(0.29 + 0.332a; - 0.331a;^), (8) 
and metal-rich (upper branch) regimes, 

12 + log{0/H\ = 12 - 2.939 - 0.2a; - 0.237a;2 - 0.305a;^ - 0.0283a;^ 
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7/(0.0047 - 0.0221a; - 0.102a;2 - O.OSlTx^ 



0.00717a;^), 



(9) 



where 



X = logR23 



log {[O I I]X3727 + [O///] AA5007, 4959) /H(3) 



(10) 



and 



y = log032 



log{[OIII]XX5007, 4959)/[0//] A3727). 



(11) 



The R23 cahbration has an estimated uncertainty of ±0.10 in log(0/H) (Kobulnicky et 
al. 1999). 

Class 3. [OII]A3727 is not detected, but [OIII]A4363 and [OIII]AA4959, 5007 are detected 
(NTT data). It is not possible to break the degeneracy using the [NII]A6584/[OII]A3727 and 
[N11]A6584/Ha ratios. 

In this case, we can derive the electron temperature but not the O /H abundance. We use 
the temperature obtained from standard nebular analysis to select the appropriate metallic- 
ity calibration with which we derive the oxygen abundance. After determining the metallicity 
regime for each galaxy, we compute the 0/H ratio using the upper branch given by Edmunds 
& Pagel (1984), 

log{0/H)^ ~ -Q.miogR^ - 3.24, (12) 

(0.6 < logR^ < 1.0) 

D _ Fyoiinw 5007,4959 
^3 = ^ 

for galaxies in the high metallicity regime, and the lower branch 



(0.4 < logR^ < 1.1) 

for low metaUicity objects. The R3 calibration has an estimated uncertainty of ±0.20 in 
log{0/H) (Edmunds & Pagel 1984). 

Two additional empirical tools were used to obtain the oxygen abundance: the P-R23 
calibration, with an estimated uncertainty of 0.10 dex (Pilyugin 2000, 2001a and 2001b), 
and the [OIII]/[NII] calibration given by Pettini and Pagel (2004), with an uncertainty of 
0.25 dex. In particular, for NTT data we use the calibration extrapolated to [OII]A3727 = 0, 
where the uncertainty increases to 0.30 dex (Pilyugin 2001a). 

To compare the derived oxygen abundances, we use as many of the methods described 



log{0/H)i ~ 1.. 



GllogRs + 6.43, 



(13) 
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Table 4. Oxygen abundance obtained using different methods. 



Method 




R Branch 


Te 


R23 


N203 


R3 


R23 Pilyugin 


Ra Pilyugin 


Adopted 


A Max 


Literature 


Fairall44 


upper 


- 


- 


- 


8.45 


- 


8.26 


8.45±0.20 


0.19 


- 


III Zwl07 


upper 


8.52ib0.12 


8.57 


8.23 


8.29 


8.31 


8.39 


8.52±0.12 


0.34 


8.30(d) 


Mrk309 


upper 




9.04 




8.98 


9.40 


10.07 


9.04±0.20 


1.09 


9.03(h) 


Mrk475 


lower 


7.94ib0.01 


8.00 


■ 


7.95 


7.81 


7.58 


7.94±0.01 


0.39 


7.93(e) 


Mrk712 


upper 




~ 




8.28 


~ 


8.37 


8.28±0.20 


0.09 


8.26(a) 


MrKlz / i 


lower 






8.22 


7.57 




8.49 


oo_Ln on 


0.92 


7.99(e) 


NGC450 


upper 


8.23±0.03 


8.61 


8.15 


8.26 


8.33 


8.33 


8.23±0.03 


0.18 




NGC4385 


upper 






8.69 


8.79 




9.61 


8.69±0.20 


0.92 


8.70(i) 


NGC4861 


lower 


7.91±0.02 


7.89 


7.92 


7.98 


7.79 


7.61 


7.91±0.02 


0.37 


7.99(e) 


NGC5430 


upper 




8.97 


8.75 


8.97 


9.08 


10.04 


8.97±0.20 


1.29 


8.81(b) 


NGC5471 


lower 


8.08±0.01 


7.99 


8.01 


7.99 


7.82 


7.61 


8.08±0.01 


0.38 


8.20(c) 


NGC6764 


upper 


8.89±0.20 


8.90 


8.78 


8.89 


8.92 


9.85 


8.89±0.20 


1.07 


8.60(b) 


NGC6764 


upper 


8.89±0.20 


8.93 


8.78 


8.89 


8.98 


9.85 


8.89±0.20 


1.07 




NGC7714A 


lower 


7.90±0.13 


7.73 


8.51 


7.05 


7.38 


6.80 


7.90±0.13 


1.71 


8.14(f) 


NGC7714B 


lower 


8.03±0.05 


8.00 


8.26 


7.65 


7.70 


7.33 


8.03±0.05 


0.93 


8.19(f) 


UM48 


upper 


8.54±0.01 


8.76 


8.44 


8.54 


8.62 


9.01 


8.54±0.01 


0.57 


8.52(g) 



(a) Contini et al. (1995), (b) Contini ct al. (1997), (c) Evans (1986), (d) Gil dc Paz et al. (2000), (e) GIT2000, (f) Gonzalez- 
Delgado et al. (1995), (g) Mascgosa ct al. (1991), (h) SGTI 2000, (i) Sugai ct al. (1992). 



above as possible (based on available spectral features) to derive 0/H values for all of the 
galaxies in our sample. Figure 5 shows the comparison between results obtained applying 
the Te-Method and the other methods used in this work. We see that in the absence of 
[OIII]A4363, use of the electron temperature Te, the R23 method (for Class 2 objects), 
and the [OIII]/[NII] calibration all provide good determinations of gas abundance. The R3 
technique was used only for MRK 712. We find that the R23 and R3 calibrations of Pilyugin 
do not perform correctly for either high or very low gas abundances, and are therefore not 
used in this work. 

Furthermore, we have compared our derived gas abundances with those in the literature 
in Figure 6. The top panel shows this comparison, while the bottom panel shows the de- 
viations as a function of 12-|-log(0/H).The median difference is 0.02 dex with a standard 
deviation of 0.17 dex. 

The above comparisons, both internal and external, demonstrate that our gas abundance 
determinations are reliable and robust. The gas abundances must be measured correctly to 
accurately constrain the massive stellar population and its relation to gas metallicity. The 
abundances obtained using the different methods and the final adopted value for the oxygen 
abundance of each galaxy are listed in Table 4. 
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5 MASSIVE STAR POPULATIONS: WR AND O STAR NUMBERS 

The number of massive stars present in a region of a Wolf-Rayet galaxy can be derived 
using the standard method developed by Conti (1991). A representative number of O and 
WR subtype stars can be obtained directly from the optical spectra (Vacca & Conti 1992, 
hereafter VC92). 

The absolute number of WR stars in the sampled galaxies was estimated using the blue 
bump (A4686) and red bump (A5808) luminosities. The blue bump is a blend of NVA4605, 
NIIIAA4634, 4640, CIIIA4650, CIVA4658 and HeIIA4686 broad WR hues (GITOO). Super- 
posed on the blue bump may be [FeIII]A4658, HeIIA4686, and [ArIV]AA4711, 4740, which 
are narrow nebular lines. The dominant contribution to the blue bump is from the broad 
HeIIA4686 line arising in WNL stars, while the red bump is mainly due to broad CIV A5808 
from WC stars (Schaerer et al. 1999, GITOO). However, some contribution from early WN 
stars (WNE) might be present in the blue bump (SV98). The detectability of the red bump 
is low compared to the blue bump. Figure 7 shows the spectral region of those galaxies in 
our sample containing either the blue or red bump. 

The dominant subtype of WR stars can be constrained by analyzing the presence of 
particular lines. N111A4640 and/or C111A4650 are observed in 9 galaxies, while NVA4604 is 
absent. The former are due only to WN stars, thus indicating the predominance of WNL 
stars (Schaerer et al. 2000). 

To derive the number of WNL stars from the blue bump, we measured the flux of the 
entire bump, and then subtracted the nebular hues and the contributions of NIIIA4640 and 
CIII/CIVA4650/A4658 when these hues were present. This procedure then retains only the 
broad emission component of HellA4686. 

The CIIIA5696 emission line is a signature of late type WC stars (WCL, GITOO), and their 
presence is expected in high metallicity regions (Maeder, 1991). When CIIIA5696 is absent, 
but CIVA5808 is present, the emission is likely due to early-type WC stars (WCE). WN 
stars cannot be responsible for the CIVA5808 emission, since the ratio HeIIA4686/CIVA5808 
observed in our spectra ranges from 0.64 to 5, much lower than the value of ~ 16 predicted 
by SV98 (Schaerer et al. 2000). 

The absolute population of WR subtype stars can be derived if the line luminosity 
of a single subtype star is known. We use the integrated broad emission line luminosity 
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of HeIIA4686 in each of our galaxies to derive the number of WNL stars, the CIVA5808 
luminosity to derive the number of WCE stars, and CIIIA5696 for WCL stars (GITOO). 
In practice, 

T obs 

J-^WNL = 7 , U4j 

-t>He//A4686 

T obs 

J^WCE = 7 , (loj 

-t^C/yA5808 

and 

T obs 

NwcL = (16) 

where L°^liixi&s%i ^civ \5808 ^ciii \5696 total luminosities observed around the 

corresponding spectral features. For the average luminosity of a single WNL star in the He 
II A4686 hue we use Lnexmse — 1-6 ± 0.8 x 10^^ erg s~^, while for a single WCE star in the 
CIVA5808 line we use Lcivxbsos = 3.0 ± 1.1 x 10^^ erg s~^, and for a single WCL star in the 
CIIIA5696 line we use Lcmx^me = 8.1 ± 2.9 x lO^^ erg s-i(SV98). 
Then, the total number of WR stars is defined by 

NwR — NwCE + NwCL + NwNL- (17) 

We note that these line luminosities are based on WR stars observed in the Milky Way, 
assuming solar metallicity . The line luminosities show significant scatter depending on the 
dominant WR subtype and the metallicity in the observed object. 

To obtain the number of O stars in each galaxy, we assume that all ionizing photons, 
Qq, are produced by O and WR stars (Conti 1991). Hence, 

Qobs ^ No7vQo7V + NwrQwR (18) 

where Norv is the number of 07V stars and Qwr and Qo7V are the number of ionizing 
photons per second produced by WR stars (all subtypes summed up) and 07V stars, respec- 
tively. Thus, the total number of ionizing photons can be obtained from the H/5 luminosity, 
L{H(3), through the relation 

g°^^ = 2.01 X 10^^Lh/3. (19) 

The number of O stars present is derived from the number of 07V stars after applying a 
correction for the presence of other O star subtypes. In this sense, VC92 and Vacca (1994) 
defined the conversion parameter as the proportion of 07V stars relative to all OV stars, 

m = No7v/Nov. (20) 
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Table 5. Star Pojjulatiori Pararnotors. 



Qo(xl052) Age (Myr) j?o(t) 



Fairall44 

III Zw 107 

Mrk309 

Mrk475 

Mrk712 

Mrkl271 

NGC450 

NGC4385 

NGC4861 

NGC5430 

NGC5471 

NGC6764 

NGC6764 

NGC7714A 

NGC7714B 

UM48 



89.76±13.46 



0.12±0.02 
2.27±0.34 
0.04±0.01 
3.52±0.53 
9.85±1.48 
3.18±0.48 
15.28±2.29 
0.03±0.00 
14.94±2.24 
12.40±1.86 
29.68±4.45 
2.79±0.42 
9.94±1.49 



3.96±0.59 
18.61±2.79 



3.0 
4.0 
6.1 
4.1 
3.5 
4.7 
2.6 
4.9 
4.0 
5.1 
3.8 
5.5 
5.5 
5.4 
3.7 
4.8 



0.50 
0.53 
0.80 
0.48 
0.86 
0.36 
1.50 
0.18 
0.53 
0.14 
0.65 
0.12 
0.12 
0.52 
0.72 
0.20 



The parameter 770 depends on the IMF for massive stars and is a function of the time 
elapsed since the beginning of the burst. 

Using the models of SV98, which give the evolution of EW(H/3) as a function of the time 
elapsed from the beginning of an instantaneous burst for different metallicities, we derive 
the starburst age t for each of our galaxies. These ages are in good agreement with the 
predicted age and duration of the WR phase estimated from instantaneous burst models. 
We use the models of SV98 to evaluate the parameter rio{t), adopting the oxygen abundance 
obtained for objects as the metallicity, the canonical slope for a Salpeter initial mass function 
r = —2.35, and a stellar mass upper hmit of 120Mq. The number of ionizing photons, the 
age of the starburst and ri{t) derived for each galaxy of our sample are given in Table 5. 

Then, using equations (18), (19) and (20), the absolute number of O stars can be derived 

as 



We adopt Qwr = Qorv = 1-0 x lO^^s-^ (Schaerer et al. 1999). 

The absolute number of O stars, Wolf-Rayet stars, ^wnl and ^wce, and the ratio 
NvTii/No are given in Table 6. 

5.1 Comparison with Evolutionary Models 

The relationship between the ^wnf^o ratio and oxygen abundance obtained for the galaxies 
in our sample is shown in Figure 8. We adopt a solar oxygen abundance of 12+ log(0/H) 
= 8.70 (Grevesse & Anders 1989, Grevesse & Sauval 1998). The predictions of Schaerer & 



No = Nov = 



Qo"^ — NwrQwr 
'noit)Qo7v 



(21) 
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Table 6. Massive Star Pojjulatiori. 



^WNL ^WCL ^WCE ^WR Nq ^Wr/I^O* ^Wc/^WN* 



Fairall44 


692±363 






692±363 


5660it623 


0.12±0.04 




III Zw 107 


3389±2211 






3389±2211 


28726±4172 


0.13±0.04 




Mrk309 


30841±16314 


7945 ±229 




38786±16315 


58498±16883 


0.66±0.24 


0.26±0.09 


Mrk475 


25±14 




7±3 


32±14 


186±29 


0.17it0.16 


0.28±0.10 


Mrk712 


661±334 






661±334 


1876±390 


0.35±0.13 




Mrkl271 


11±7 






llit7 


77±22 


0.14±0.05 




NGC450 


322±215 




129±51 


451±221 


2044±147 


0.22±0.08 


0.40±0.14 


NGC4385 


3037±1534 


1208±960 


972±634 


5217±2288 


25740±7610 


0.20±0.07 


0.72±0.25 


NGC4861 


522±288 


65±3 


58±35 


586±290 


4375±489 


0.15±0.05 


0.24±0.08 


NGC5430 


4761±2520 




999±435 


5761±2558 


66149±17776 


0.09±0.03 


0.21±0.08 


NGC5471 


2±1 






2±1 


44±2 


0.05±0.02 




NGC6764 


2812±1687 


859±42 


1750±906 


5421±2231 


79325±14771 


0.07±0.02 


0.93±0.33 


NGC6764 


2380±1457 


735±49 


1606±817 


4721±1943 


56180±11250 


0.08±0.03 


0.98±0.35 


NGC7714A 


1920±1288 


869±18 


1596±883 


5295±2277 


44282±2655 


0.12±0.04 


0.87±0.31 


NGC7714B 


451±255 






451±255 


3276±357 


0.14±0.04 




UM48 


1654±1279 




353±234 


2007±1300 


38967±6385 


0.05±0.02 


0.21±0.08 



Vacca (1998) models for instantaneous bursts (solid line) with IMF slopes of -1, -2 and -2.35 
(Salpeter), and of Starburst 99 (Leitherer et al. 1999) for extended bursts (duration of 2-4 
Myrs) with IMF slope = -2.35 and mass hmit of 100 Mq (dashed line) are plotted. 

Evolutionary models predict that for a given metaUicity, the ratio between WR and O 
stars varies strongly with the age of the burst, and the duration of the WR stage in the 
starburst also increases with increasing metaUicity. (Maeder & Meynet 1994; Meynet 1995; 
SV98). Maeder (1991) interpreted this behavior as the result of increased stellar mass loss 
at higher metallicities. The increased stellar mass loss reduces the mass hmit for forming 
WR stars in metal rich galaxies (Maeder 1991). 

In the low metaUicity case, the results can be explained by an instantaneous burst with a 
Salpeter IMF slope (F = -2 to -2.35). In the high metaUicity regime, the results deviate from 
the expected behavior based on models with an instantaneous burst for starburst galaxies 
(SV98). A steeper IMF slope is required than in the low metaUicity case, or we must invoke 
an extended burst to explain the results. SGITOO interpret this behavior for their sample of 
high metaUicity galaxies as implying an extended burst duration of ~4-10 Myrs. Support for 
this conclusion comes from the observed WR population relation (WC/WN) with metaUicity 
and the red supergiant star features observed in their objects (GITOO, SGITOO). NGC4385, 
NGC 5430, NGC 6764 and UM 48, with 12 + log{0/H) > 8.4, are candidates for this type 
of object. 

In Figure 9 we plot the ratio ^wcf^WN (where WC = WCL -|- WCE) versus oxygen 
abundance compared to the predictions of SV98 models for an instantaneous burst, as well 
as extended bursts of 2 and 4 Myrs (Pindao et al. 2002). The dot-dashed hue (bottom right) 
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shows the observed trend of WC/WN with metaUicity in Local Group galaxies derived 
empirically from observations by Massey & Johnson (1998). The ^wcf^WN ratio is lower 
than the instantaneous burst model predictions and, for galaxies with metaUicity higher 
than 12 + log{0/H) > 8.4, the ratio is closest to the values reported for galaxies in the 
Local Group corresponding to a constant star formation regime (Massey & Johnson 1998). 

These low ^wcf^WN ratios might be partially explained by the assumption that the 
predominant contribution to the CIVA5808 line luminosity at high metallicities is from 
WGE stars. If instead we assume that the main contribution to this line luminosity comes 
from WGL stars, the lower luminosity of WGL stars can increase the ratio of WG to WN 
stars by a factor of ~ 3 — 4 (GITOO). Schaerer & Vacca 1998 predict a luminosity ratio 
GIIIA5696/GIVA5808 for a WGL (WG7) star of ~ 0.5, but galaxies in our sample where 
WG emission lines are measured show lower GIIIA5696/GIVA5808 ratios. 

We contend that the low Nwcf^WN ratios observed in high metaUicity galaxies are due 
to the nature of the star formation bursts, with different durations for the WN and WG 
stages (GITOO). These galaxies have low H/3 equivalent widths leading to ages greater than 
5.3 Myr (except UM48 with age ~ 4.8 Myr). They are therefore in the late stages of their 
WR episodes (GITOO). At these later times, WN stars are still present, while the number 
of WG star drops to zero (SV98). This supports an extended burst to explain the N^yjv/No 
ratio seen in these galaxies. 

The presence of late type stellar features provides additional evidence supporting an 
extended burst. These absorption features are commonly observed in integrated spectra of 
stellar clusters, indicating the presence of red giants and supergiants (Bica & AUoin 1986). 
The presence of late-type stars further supports the idea of an extended burst with an age 
<7 Myr. In Figure 10 we present the spectra of these four high-metallicity galaxies which 
show apparent TiO bands (~ A6250) and and a blend of FeI+BaII+GaIA6495, which are 
spectral features of characteristic of latc-typc stars. 

Alternatively, Pindao et al. (2002) suggest that corrections to the evolutionary models 
can reproduce the low observed ^wr/^o ratios using an instantaneous burst. They invoke 
the uncertainties in synthesis models of WR bumps to explain their results, since the WR 
line luminosities show a large scatter in the WNL calibration sample of SV98 with Galactic 
and LMG objects. A dependence between the WR luminosity and the bolometric stellar 
luminosity may cause this discrepancy. The WR line luminosity L4686 is observed to increase 
with increasing bolometric stellar luminosity. They suggest a more precise analysis, splitting 
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the bolometric stellar luminosity in two domains {log[L/ Lq) < 6 and log{L / Lq) > 6). The 
different average line luminosities for each domain are L4686 = 5.6 x 10^^ erg s~^ for log{L) < 6 
and L4686 = 3.1 X 10^^ erg for log{L) > 6. Using these values in synthesis models leads 
to a reduction in the predicted number of WNL stars. Note that the hue luminosity for a 
WNL star with L > 6Lq is almost twice as large as that adopted by SV98 for their Galactic 
cahbration sample. Using these values in the synthesis models, their observed WR features 
can be reproduced with an instantaneous burst and a standard Salpctcr IMF (Pindao et al. 
2002). The recalibrated (lower) luminosity of WN stars used in models leads to a significant 
reduction in the WR bump intensity in the bursts with ages > 4-5 Myrs. 

Only the youngest bursts with very high EW{H(3) are dominated by very luminous 
WNL stars. For objects with low N^r/Nq ratios, if we apply the suggested correction to 
the L4686 luminosities in the actual observations (instead of the models), and recalculate 
the number of WN stars at the highest metallicities (dominated by low luminosity WNL 
stars), the lower luminosity of WNL stars can increase the fraction of WN stars by a factor 
of ~ 2 — 2.3 (GITOO). With this change in assumed luminosity, we are unable to reproduce 
the standard instantaneous burst model predictions for the Nwr/Nq ratio in the highest 
metaUicity galaxies. This analysis confirms that the low values of ^wr/^o for objects with 
12 + log{0/H) > 8.4, obtained for NGC4385, NGC 5430, NGC 6764 and UM 48, require an 
extended burst with a Salpeter IMF slope. 

6 CONCLUSION 

In this work, we present a spectroscopic study of 14 Wolf-Rayet galaxies from the sample of 
SCP98, as well as NGC450 for which the WR features are newly detected. Our goals were to 
search for and confirm the presence of WN and WG stars, and to compare the results with 
predictions from evolutionary synthesis models (SV98). We tested the agreement of these 
models with observations for a large range of metaUicities, spanning 7.90 < 12 -|- log(0/H) 
< 9.04. 

Our main results can be summarized as follows: 

1. The broad WR emission in the blue region of the spectrum, the blend of NIIIA4640, 
GIIIA4650, GIVA4658, and HeIIA4686 emission lines, is present in all fourteen galaxies. The 
WR population in these galaxies is dominated by late WN stars. However, the red bump 
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produced by the emission of broad C IVA5808 from early WC stars is detected in only nine 
galaxies. 

2. The weak, broad WR emission line CIIIA5696 is detected in six galaxies, which suggests 
the presence of late WCL stars in these objects. This line is expected in high-metaUicity envi- 
ronments (GIT2000). A good example of this is Mrk 309, with a metaUicity of 12+log(0/H) 
= 9.04, where this line is strong. 

3. We found good agreement when comparing the relative numbers of WR and O stars 
{Nwr/Nq) obtained from observations and those predicted by the evolutionary synthesis 
models of SV98 for low metaUicity galaxies. The ratio Nwr/Nq in these galaxies can be 
explained by a burst of star formation. The Nwr/Nq value and the observed equivalent 
widths of the blue and red bumps also compare favorably with SV98 predictions. We found 
that it is necessary to invoke an IMF slope between — 2 < F < — 2.35 and an instantaneous 
star formation event to explain the observed Nwr/Nq ratios in low metaUicity regimes. 

4. For NGC4385, NGC 5430, NGC 6764 and UM 48, the Nwr/Nq ratios are lower than 
the predictions of models with an instantaneous burst for metallicities of 12 + log{0/H) > 
8.4. The existence of an extended burst is supported by the presence of TiO bands (~ A6250) 
in these objects spectra. The presence of this older stellar population is indicative of the 
necessary time elapsed since the burst for stars to evolve through this phase. The Nwn/Nwc 
ratios combined with their high metaUicity suggest these objects are in the late stages of 
their WR episodes (> 5.3Myr), when WN stars arc present while the number of WC stars 
drops to zero (SV98). The observations can be better represented using extended starbursts 
with durations of 2-4 Myr and Salpeter IMF slope. 

5. Massive stellar evolution models predict that the relative number of WR stars increases 
when metaUicity increases (Maeder & Meynet 1994, SV98). Our results do not confirm this 
trend. The partial disagreement between our results and the models may arise from large 
uncertainties in the luminosity of a single WR star and the uncertainties associated with the 
best choice of the dominant contribution of WR star subtype in the high metaUicity regime. 

The massive stellar populations in Wolf-Rayet galaxies can be better understood if stud- 
ies like this one are combined with other stellar population investigations. A detailed spec- 
tral study in the infrared of molecular bands and low ionization absorption lines can provide 
information about older stellar populations in a starburst (Origlia et al. 1999, GITOO). Spec- 
troscopy in the UV (912 — I8OOA) can provide information about the young star population 
in galaxies. This UV spectral range contains resonant spectral lines of OVIA1035, SiVA1400 
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and/or CIVA1550 that are spectral signatures of young massive stars (Leitherer et al. 2002). 
The combination of the IR and UV yields information on the current stellar populations, 
and constrains the upper and lower limits of time elapsed since the burst (Leitherer et al. 
2002, Olivia et al. 1999). Additionally, analysis of the ct-element abundances can provide 
additional constraints on the age of the burst and consequently the number of massive stars 
(Lanfranchi, G. A. & Priaga 2003). 
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APPENDIX A: REMARKS ON INDIVIDUAL OBJECTS 

In this section we present a brief description of some specific properties of each galaxy in 
our sample. By doing so we can address specific issues related to these objects and compare 
with results from other authors, and obtain important information for the analysis of the 
population whole. 

Fairall 44: Kovo & Contini (1999) reported the blue bump in their systematic search for 
Wolf-Rayet stars in young starburst galaxies. We report the presence of WN stars in this 
galaxy. We find that the WR/0 ratio is 0.24. 

Ill Zw 107: Kunth & Jobert (1985) reported a moderately strong emission band at 
4686 A due to WR stars in III Zw 107 S. We report the presence of NIIIA4640 and HeIIA4886 
broad emission lines in this galaxy. The WR/0 ratio is 0.12. 

Mrk 309: This bright UV continuum galaxy has broad emission features at NIIIA4640 
and HeIIA4686 from WR stars noted by OC82. They found that the number of WR stars 
is comparable to the number of O stars in this galaxy and the nuclear region of NGC6764. 
We note the presence of NIIIA4640 and HeIIA4686 in the blue bump and CIIIA5696 in the 
red bump in our spectrum of this galaxy. We find that the number of WR stars is similar 
to that of O type stars in the nuclear region of Mrk 309. Our results are in good agreement 
with SGITOO. 

Mrk 475: Broad HeIIA4686 and NIIIA4640 emission hues were first noted by Conti (1991). 
Strong blue and red bumps were detected by ITL94. GITOO identified features of SiIIIA4565, 
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HeIIA4686, HeI,NIIA5047, and CIVA5808, where the blue bump is strongly contaminated by 
nebular emission. We confirm the results of GITOO, and find a WR/0 ratio of 0.17. 

Mrk 712: Contini et al. (1995) reported the discovery of emission from WR stars in the 
giant HII region 4.5" south of the nucleus. They estimated a WN/0 ratio of 0.2 from the 
HeIIA4686 luminosity. The [ArV] emission hue in their spectrum showed that the HII region 
is strongly ionized by hot WR stars. The WN/WO ratio indicated a very young starburst 
episode and a flat initial mass function with slope F between -1 and -2. We use the SV98 
evolution models to compare with our observational results (§6.2). 

Mrk 1271: Izotov & Thuan 1998 detected a broad blue bump. GITOO identified the 
NIIIA4512, SiIIIA4565, NIIA4620 and HeIIA4686 emission hues, and possibly CIVA5808 in 
their high quality spectrum. They reported the blue bump as being strongly contaminated 
by nebular emission. We confirm only the presence of the HeIIA4686 broad line in the nuclear 
region of this galaxy. The CIVA5808 broad line is not detected in our spectrum of Mrk 1271. 
The WR/0 ratio is 0.14. 

NGC 450: There is no mention of WR detection in this galaxy in the literature. We report 
a ratio of WR/O=0.22. The models of SV98 for an instantaneous star formation episode 
predict a flat IMF slope (F = — 1 to —2) to explain the presence of WR stars in this low 
metallicity galaxy. 

NGC 4385: HeIIA4686 and NIIIA4640 Wolf-Rayet features are present in this starburst 
galaxy, according to Durret & Tarrab (1988). An optical spectrum by Salzer (1990) shows 
essentially the same features. Conti (1991) identified a narrow emission feature near 4660A as 
[Felll], while Salzer et al. (1989) identifies it as CIV. We detect NIIIA4640 and HeIIA4686 
in the blue bump and CIIIA5696 and CIVA5808 in the red bump. We find a high number 
of WR stars compared to O stars. The WR/0 ratio is 0.13 for the nuclear region of NGC 
4385. 

NGC 4861: Dincrstcin & Shields (1986) and Izotov, Thuan & Lipovetsky (1997) detected 
the blue and red bumps. GITOO identified NIIIA4512, SiIIIA4565, 7Vl^A4619, HeIIA4686 and 
CIVA5808 broad emission lines, and a blue bump strongly contaminated by nebular emission. 
Even with their high-quahty spectrum, the NIIIA4640 and CIVA4658 hues appear blended. 
In our spectrum of the central region of this galaxy, we can resolve the lines of the blue 
bump. We report the presence of CIIIA5696 and CIVA5808 broad emission hues due to WC 
stars. 

NGC 5430: A strong emission feature near 4650A due to WR stars in the spectrum of 
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a bright region SE of the galactic center was noted by Kell (1982). He identified NIIIA4640 
and HeIIA4686 emission lines as coming from WN stars (Kell 1987). No strong emission 
features are seen in the lUE spectrum of this source. A relatively older stellar population 
appears in the center of NGC 5430. According to Kell (1987), the knot 20" southeast of the 
nucleus might be a separate galaxy interacting with NGC 5430. For this knot we measure 
the NIIIA4640 and HeIIA4686 blue bump broad components, as well as the CIIIA5696 and 
CIVA5808 broad emission lines. A high WR/0 ratio of 0.09 is found for NGC 5430. 

NGC 5471: This is a massive giant HII region located in MlOl. Mass- Hesse & Kunth 
(1991) found WR-HeIIA4686/H/3 = 0.02 and a 3.5 Myr burst. NGC 5471 seems to be dom- 
inated by a well-defined burst of star formation. A broad HeH emission fine (FWHM~2000 
km s~^) was detected in this region by Castaneda et al. (1990) and was confirmed by Mass- 
Hesse et al. (1991). We find a low ratio WR/O=0.04 for this galaxy. 

NGC 6764: OC82 noted the presence of broad NIIIA4640 and HeIIA4686 emission fine 
features from WR stars. They attributed the 4660A emission line to GUI and not to [Felll] 
forbidden emission. The similarity of the overall emission line spectrum of NGC 6764 to that 
of dwarf galaxies and to certain giant HH regions attracted the attention of these authors. 
We observed this object at two position angles, 90° and 46°. We find WN and WC stars in 
the nuclear region of this galaxy. For both position angles the number of Wolf Rayet stars 
is comparable to that of O type stars. 

NGC 7714: Weedman et al. (1981) called this object a "prototype starburst" galaxy, van 
Breugel et al. (1985) noted the spectral similarity of NGC7714 to Minkowski's object and 
extragalactic HH regions. They reported this galaxy as having weak WR emission features 
near HenA4686. Conti (1991) called attention to the importance of NGC 7714 for under- 
standing the relation between starbursts and the presence of WR stars in this kind of galaxy. 
We found WR stars in two regions of NGC 7714. In the nuclear region we find WR/0 = 
0.08, while in the secondary HII region, WR/O=0.14. In the nuclear regions we were able 
to measure the CIIIA5696 and CIVA5808 broad emission fines due to WC stars. 

UM48: A systematic search for Wolf-Rayet features in this galaxy was done by Masegosa 
et al. (1991). The authors reported the presence of a WR blue bump. They assumed that 
the global WR detection rate depends on the metaUicity and they analyzed the SN He 
contribution to the blue bump. They concluded that a large blue bump luminosity should 
be expected in regions with SN contamination. We did not detect the NIIIA4640 broad 
emission line in this galaxy. We measure the WC broad emission line CIVA5808. 
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(g) (h) 

Figure 1. DSS images of galaxies with slit positions, (a) Fairall 44, (b) III Zw 107, (c) Mrk 309, (d)Mrk 475, (e) Mrk 712, (f) 
Mrk 1271, (g) NGC 450, (h) NGC 4385. The horizontal bar at the lower left indicates 0.5 arcminutes. 
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(o) 

Figure 1. Continued. DSS images of galaxies with slit positions, (i) NGC 4861, (j) NGC 5430, (k) NGC 5471, (1) NGC 6764 
P.A. 67, (m) NGC 6764 P.A. 90, (n) NGC 7714 and (o) UM48. The horizontal bar at the lower left indicates 0.5 arcminutes. 
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Figure 2. Top: Comparison of the values of E(B-V) adopted in this work and those obtained by other authors. The dashed 
Une shows a relation of unity. Bottom: The difference between the values adopted in this work and those obtained by other 
authors. 
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Figure 3. The diagnostic diagram of [OIII]A5007/H/3 plotted against [OI]A6300/Ha.The solid line shows the border between 
starburst and AGNs ionization mechanisms. The dashed lines show the model uncertainty of 0.1 dex. 
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Figure 4. The diagnostic diagram of [OIII]A5007/H/3 plotted against [iV/i']A6584/HQ:.The solid line shows the border between 
starburst and AGNs ionization mechanisms. The dashed lines show the model uncertainty of 0.1 dex. 
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Figure 5. Compajrison between the values of 12 + log(0/H) obtained with the Te method and other methods used in this 
work. The median offset and c for each relation are noted in each panel. The R23 and R3 calibrations of Pilyugin are shown 
separately on the left. 
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Figure 6. Top: Compajrison of the values of 12 + log(0/H) adopted in this work and those obtained by other authors. The 
dashed Une shows a relation of unity. Bottom: The difference between the values obtained by other authors and those adopted 
in this work. 
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Figure 7. Enlargement of blue and red bump spectral regions of galaxies where the broad lines HeIIX4686 and CIVA5808 are 
visible. 
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Figure 7. Continued. 



© 2004 RAS, MNRAS 000, 1-25 



Massive Star Populations in Wolf-Ray et Galaxies 37 




Figure 7. Continued. 
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Figure 8. Nw r/ versus 12 + logiO/ H) for the galaxies in our sample. Predictions of Schaerer &; Vacca (1998) models for 
instantaneous bursts (solid line), and of Starburst 99 (Lcithcrer ct al. 1999) for extended burts of star formation of 2 and 4 
Myrs. (dashed line) are overplotted. The lines are labeled with the IMF slopes. 
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Figure 9. Nwc / N versus 12 + logiO/ H) for the high mctaUicity galaxies in our sample. Predictions of Schaerer & Vacca 
(1998) models for maximum value for instantaneous bursts (solid line), and SGITOO extended bursts of 2 - 4 Myrs (dashed 
lines) are overplotted. The dot-dashed line (bottom right) shows the observed trend of WC/WN with metallicity in Local 
Group galaxies derived empirically from observations by Massey & Johnson (1998). 
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Figure 10. Spectral region showing TiO bands (~ A6250) and narrow Fel+Ball+Cal A6495 for the high-metallicity galaxies 
NGC 4385, NGC 5430, NGC 6764 (P.A. 67 and 90) NGC 5471 and UM 48. 
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